The Design of Radar Signals Having 

Both High Range Resolution and 

High Velocity Resolution 
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(Manuscript received April 5, 1960) 

Radar systems that yield simultaneous inform,ati(m regarding the range 

and velocity of a target would be useful in certain applications. A discussion 
is presented of soine fundamental limitations on simultaneous range-velocity 
determination in terms of ambiguity diagrams, with the aid of a quantum 
mechanical analog to iheWigner distribidion function. A sequence of sig- 
nals is found m which the signals yield both range and velocity information 
with increasing accuracy. However, this desired property is not accompanied 
by a waveform suitable for maximum operating efficiency. 

I. INTRODUCTION 

The vciy large power requirements of present-day radars clearly de- 
mand full utilization of the potential of existing equipment. Unfortu- 
nately, an efficient equipment use does not always lead directly to a sig- 
nal possessing desired high resolution characteristics. The "Chirp" 
scheme, which employs linear frequency modulation and is discussed in 
detail in the accompanying paper/ represents one method by which the 
desirable high resolution properties may be secured with an optimum 
equipment utilization. However, in Section 3.5 of the accompanying 
paper it was pointed out that the Chirp scheme, like ordinary radar 
techniques, possesses an inherent ambiguity in a simultaneous deter- 
]nination of both the range and velocity of a moving target. Ambiguity 
diagrams," which are representations of joint response functions in both 
time and Doppler frequeiicy, provide a pictorial representation of the 
uncertainty in determuiing range (time) and velocity (Doppler fre- 
quency). Fig. 1, (which is a duplicate of Fig. 38 of the accompanying 
paper) schematically illustrates the ambiguity diagrams corresponding 
to three different transmitted signals: (a) a long constant-frequency sig- 
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Fig. 1 — (a) Schematic, ambiguity diagram for a long constant-frequency radar 
signal; the more shaded portions denote regions of greater signal response with 
time, t, and for variously moving targets with Doppler frequency', v. (b) Hchematic 
ambiguity diagram for a short constant-frequency radar signal; improved time 
resolution is secured at the expense of some discrimination against moving targets, 
(c) Schematic ambiguity diagram for a typical Chirp signal ; an approaching target 
with a positive Doppler shift has a response attentuated in amplitude and shifted 
to an earlier time. 
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nal, (b) a very short constant-frequency signal and (c) a Chirp signal 
with linear frequency modulation. 

For certain applications it may be desirable to transmit a signal whose 
ambiguity function is highly peaked only about the point ( ^ c ^ 0. 
If such a signal exists, it would permit a high resolution in both range 
and velocity simultaneously. In Section II of this paper some proper- 
ties of ambiguity functions are discussed with the aid of quantum me- 
chanical analogs, and one interesting general property is discussed in 
detail: If a requirement of rotational invariance is imposed on ambiguity 
functions (a concept made precise in Section II) a sequence of signals 
is obtained for which the peaking of the ambiguity function is steadily 
increased. While the signals in this class achieve the desired behavior 
in the ambiguity function, the transmitter waveforms are not uniform 
in amplitude and, therefore, cannot be transmitted with the same effi- 
ciency as the uniform amplitude Chirp signal. The notation used in 
the following analy.sis is the .same as that in the accompanying paper. 

IL SOME GENERAL PROPERTIES OF AMBIGUITY DIAGRAMS 

For convenience in the analysis, only a single target, which gives rise 
to a Doppler shift v = w/rf/c is studied, where v is the target velocity, 
fiiP is the radio-frequency carrier and c is the velocity of light. If the 
spectrum of an arbitrary complex transmitted signal is denoted by c{/) 
then the received signal becomes e(/ — v). Only receivmg networks that 
optimize the signal-to-noise ratio will be discussed. For an arbitrary 
complex signal e(0, this optimum filter is the so-called matched filter for 
that signal. If e(/) denotes the Fourier transform of the signal e(0, then 
the matched filter 7m(/) may be defined^ as i*(/). With a matched filter 
in the receiver, the response spectrum is determined by i*(/)e(/ — j^). 
The Fourier transform determines the time response: 

J — CO 

which is written as a two-parameter function. This is the expression given 
by Woodward,- and can be used perfectly well to specify the ambiguity 
diagrams. However, for the analogy to be drawn to single-particle non- 
relativistic quantum mechanics, it is convenient to develop an alter- 
nate expression for the ambiguity function amplitude. 

Suppose that the transmitter spectrum, e(/), is concentrated about 
some high rf frequency, so that 

Hf) =&{f-M, (2a) 
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where £(/) is a spectrum centered about zero frequency. It follows from 
(2a) that 

e(t) - 8(0e'"'^^^'. (2b) 

Therefore, S(0 represents the generalized envelope of the transmitted 
signal. If (2a) is substituted into (1), one readily finds that 

x(i,.) = e^"<'»-+-"" £ S* (/ + ^-) g (/ - 5) e"" df, (3a) 

or, when expressed in terms of &{t), 

Let the carrier in the presence of a Doppler shift be taken as /rf + i'/2. 
With this choice of carrier, 

,m = £ 8* (. - 6 (. + e- dr (4) 

is defined to be the complex two-parameter response envelope. For con- 
venience, the Doppler angular frequency 9, = 2-kv has been introduced 
in (4). It is the absolute value of p{t,9), which represents the detected 
output. The remaining analysis uses the definition made in (4); it is also 
the form studied by Gabor^ and Ville.^ 

No loss of generality is made by assuming the transmitted signal car- 
ries unit energy : 

r \M)\^dt = f \z{t)\^ dt = r \s(f)fdf ^ 1. (5) 

With this normalization, it follows from (4) that 

p(0,0) = [ iS(T)f rfr - 1. (6) 

This result represents a constraint that requires a portion of the "am- 
biguity" to remain in the vicinity of the point t = U = 0. By a, straight- 
forward integration using (4) and (5), the following result may be 
established : 

^ [ f {p(t,Q)\'dtdn = 1. (7) 

Zt •'-00 •'-=0 

[This result applies equally well to a similar integration involving 
xif;'')-] The interpretation given to (7) is as follows: There is, quali- 
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tatively speaking, a total ambiguity equal to "one", which must be 
distributed throughout the t,U plane. The approximately equal amounts 
of shaded area in Fig. 1 are intended to be a reflection of this constraint. 
This result seems to preclude any simple solution to achieve simultane- 
ously an arbitrarily high resolution in both range and velocity. 

Another observation worth making about the three diagrams in Fig. 1 
is that it appears they could be obtained from one another by a "rota- 
tion" about the origin of the t,Q plane. It will be demonstrated below that 
this is a general result. That is, for every ambiguity diagram defined by 
p(t,U) there exists another signal that generates another diagram differ- 
ing from the first only by a rotation about t = U = 0. The use of the 
word "rotation" needs some explanation here. Actually, the rotation can 
proceed along arbitrary elliptical contours. However, for convenience, 
it is assumed that the scales of time and Doppler frequency have been 
chosen so that the rotation is nearly circular. 

The rotation theorem, as well as others like it, is most readUy estab- 
lished by using methods common in the study of quantum mechanics. 
One of the reasons for choosing the form of (4) in preference to that of 
(1) was to facilitate the comparison with a quantity studied previously 
in quantum mechanics. Only a brief sketch is presented here of the com- 
parison used to obtain most of the remaining results of this section. 

The quantity of interest in quantum mechanics is the so-called Wig- 
ner^ distribution function defined by 

Here, h = 2Trfi is Planck's constant and )p(y) represents the Sclu'odinger 
wave function for a single particle in one dimension. Of more direct 
interest is the characteristic function for this distribution, i.e. its double 
Fomier transformation ; 

= /,*(, -|)-V(. + I) ^.. 

If this result is compared with (4), one notes an exact functional equiva- 
lence, with t playing the role of y and U playing the role of K/fi. Equation 
(9) is perhaps easier than (8) to study, and has been treated in some de- 
tail by Bass,^ MoyaF and Groenewold;^ further references may be found 
in these works. The relative simplicity of (9) stems from its simple repre- 
sentation in the quantum mechanical Hilbert space; 
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p{y,K) = (%e''''^'-'''"'^) . (10) 

Here * denotes an arbitrary normalized vector and q and p are the usual 
operators satisfying qp — pq = ifi. Numerous additional relations may 
be found, especially in the work of Moyal. 

The necessary signal transformation to generate a rotated ambiguity 
diagram will not be derived; the result wiU be stated without proof. 
Suppose that S(/) is an envelope which generates a particular ambiguity 
function p(t,il) according to (4). The modified envelope, &e(t), must be 
chosen so that its ambiguity function is defined by 

pe{t,v) = p (f cos 5 + - sin 9, J- cos d— yt sin e\ . (11) 

Here, 7 represents the ratio of the "angular frequency" axis to the 
"time" axis of the ellipse of rotation and 6 represents the "angle" of 
rotation. When j 7^ I, rotation is along nearly circular orbits. The solu- 
tion for Zeil) is given by 

Ut) = f_ U(t,T-m(T)dT, (12) 

where the kernel U{t,T\B) is expressed as follows; 

^('■^^^^ - 1/2-rik-. ^"P f^ [(^) ™* « - (r CSC .]} . (13) 

Equation (13) follows from the properties of the well-known harmonic 
oscillator problem, and is just one of the relations readily obtained with 
the techniques of quantum mechanics. 

The particular relations expressed in (12) and (13) are perhaps not 
as interesting as an additional result that can be derived from them. 
In the discussion of the ambiguities present in the interpretation of the 
diagrams in Fig. 1 it seemed that the resolution was relatively good in 
one direction but was poor in the perpendicular direction. It is just such 
a variation with direction that is not wanted in attempting to achieve 
a signal whose ambiguity function is concentrated about the point t = 
S2 = 0. For clarity, assume again that the t and 9. axes have already been 
scaled for the specific system under study. Then, to secure (as nearly 
as possible) an ambiguity diagram concentrated about the origin, it 
appears desirable to seek a rotational invariance of the ambiguity function; 
i.e., pe{t,ii) = p{fM)- It is a simple matter to demonstrate that peii,^) = 
p{t,Q.) if and only if 

£,(() = f U{t,T;6)&{r) dr - e-'^'sd), (14) 

•'—03 
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where X represents a suitable constant. There are a clenumerablc number 
of sokitions to (14), a typical one being denoted by &„{t), where X„ = 
7t -\- ^. These solutions, familiar as the eigenf unctions of the harn-ioni(! 
oscillator, are defined by the following relation: 



£,,(0 ^ 



7' 



V7ri/i!2" 



H.Wyt)c 



~A.-T'^/2 



(15) 



where Hn{z) represents the nth Hermite polynomial defined, for example, 
by 



H„(.) = (-i)v^|;^.- 



(16) 



Figs. 2(a), 3(a), 4(a) and 5(a) show four curves of &„{l) for n = 0, 4, 10 
and 2(i, respectively; each curve is symmetric about i = and is only 
shown for positive values of /. Note also a change in scale in the differ- 
ent figures. In each example 7 was chosen equal to one. 

It is a well-known fact that the functions 8„(/.) defined by (15) have 
spectra, £„(/), which are specified by the very same functional form.f 
Consequently, Figs. 2(a) through 5(a) also represent the frequency spec- 
tra for the same n values. One can see from these four figures that, as n 
increases, both the pulse width and bandwidth increase. The pulse-width- 
bandwidth product increases as {n + \y as n is increased. One may 
hope, therefore, that the ambiguity function becomes more nearly con- 
centrated about the origin as n increases. The normalization condition 
(7) must still be satisfied, but, for large n, the remaining "ambiguity" 
is thinly spread out over a large region of the i,S2 plane. 
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Fig. 2 — (a) Envelope of the transmitted signal when the parameter n = 0; in 
this ciise the curve is gaiis.sian. (b) Envelope of the "rotationally invariant" 
matched-Iilter response when h = 0; here the curve is also gaussian. 

t See Pair 702 in the tables of Campbell and Foster." 
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Fig. 3 ■ — Envelope of (a) transmitted signal and (b) matched-filter response 
when n = 4. 



This interpretation may be verified in detail by evaluating tlie am- 
biguity function p„(i,i2) when 6(0 is given by S„(0- In this calculation 
it is advantageous to employ the result that Pn{t,^) will, by construction, 
be rotationally symmetric. Hence, it suffices to study the cross section 
of the ambiguity diagram when S2 = 0; the complete contour may be 
reconstructed by rotation about the point i = S2 = 0. It is an interesting 
fact that the kernel U{t,7\6) can be rapidly converted into a generating 
function for the various response signals p„(^,0). The details of such a 
calculation are essentially present in the work of Moyal.' Either from 
the generating function or by a direct use of the definition for p(/,f2) 
one finds the zero Doppler response to be 
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Fig. 4 — Envelope of (a) transmitted signal and (b) matched -filter response 
when n = 10. 
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where Z/„(2) represents the nth Laguerre polynomial defined by 

(72" 



(17) 



(18) 



The output signal /j„(/,0) is illustrated in Figs. '2(b), 3(b), 4(b) and 5(b) 
for ft = 0, 4, 10 and 2C, respectively. Since these functions are clearly 
symmetric only positive / values have been included (note the change 
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Fig. 5 ■ — Envelope of (a) transmitted signal and (b) matched-filter response 
when n — 26. 
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of scale in the different figures). Agaiu, these curves were obtained 
assuming y = 1, In order to illustrate the rotationally invariant charac- 
ter of these signals, Fig. G shows a cut-away of a three-dimensional pro- 
jected view of the complete ambiguity amplitude p(/,S2) corresponding 
to the case n = 10. 

Fig. 6 illustrates the tendency of the ambiguity f miction to peak about 
the origin of the t,^ plane. Fmihermore, it demonstrates that a narrow 
peak is achieved only when the remaining "ambiguity" is spread out 
over a relatively' large domain of the t,il plane. 

If input and response envelopes (6„(/) and p„(/,ii)) with higher values 
of n are considered, the ambiguity function becomes more sharply peaked 
and more spread out. A study of the behavior of the zero Doppler re- 
sponse for large \'alues of n shows that 
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Fig, 6 — Cut-away of three-dimensional projected view of complete ambiguity 
amplitude p((,S2) when n = 10. 

p.C^O)- Jo(-s/2^0, (19) 

where ./o is the zeroth-order Bessel function. This asymptotic form is 
valid when / is greater than->/n but not so large that ( becomes com- 
parable with n itself. Equation (19) illustrates that, as n increases, the 
form of the response does not markedly change, only the scale changes. 
While the "rotationally invariant" response signals discussed here 
possess a narrow 3-db pulse width, they are endowed with long oscilla- 
tory tails, which fall off approximately as t~'. This rate of fall-off is to 
be compared with the rate of fall-off of the Chirp signal — namely, as 
r\ which follows from (19) of the accompanying paper. Therefore, in 
securing an ambiguity function highly peaked about the origin a penalty 
must be paid in the rate of fall-off of the time "side lobes." Nevertheless, 
as the parameter n increases there is a change of scale and the side lobes 
are pulled inward. Consequently, at least in principle, a signal S„(0 
could be chosen that mainttiined side-lobe levels below an arbitrary 
level outside a prescribed ellipse in the t,U plane. It is interesting to 
speculate whether there exists some signal envelope £(() that could be 
efficiently transmitted in addition to having a highly peaked ambiguity 
function. Unfortunately, this question remains unanswered. 
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